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Abstract 

The SABR model is shortly presented and the volatility swap explained. The fair value for a 
volatility swap is then computed using the usual theory in financial mathematics. An analytical 
solution using confluent hypergeometric functions is found. The solution is then verified using Rama 
Cont's functional calculus. 



1 Introduction 

Mathematical models for the computation of financial derivatives have gained some importance in the 
financial industry. A breakthrough vi^as realized vi^ith the Black & Scholes model. Though widely applied, 
this model soon began to show its limits and shortcomings. Among other deficiencies, the hypothesis of 
constant volatility for the underlying appeared as a much to strong assumption. It is usually now dealt 
with by quoting a so-called implied volatility, which, when inserted in the Black & Scholes equation for 
the option pricing, gives the traded prices for the latter. 



1.1 The SABR Model 

As another Ansatz, models with time-dependent volatilities or even stochastic volatilities where intro- 
duced. A famous example is the Heston-model, where the variance of the underlying is itself stochastic. 
Another widely used model is the SABR-model, where this time the volatility is stochastic. Under the 
forward measure, the forward prices Ft for an underlying is governed by the stochastic process 

dFt=atFtdWt, 

dat = aatdZt, (1) 

where a is some positive constant and where Wt and Zt are Wiener processes with a correlation coefficient 
— 1 < p < 1 [Ij. In general, the first stochastic difl:erential equation reads dFt — atFfdWt, with < /3 < 1, 
but we are going to work with /3 = 1 throughout our computations. This will have no infiuence in the 
following, since in our definition of the volatility swap, only the stochastic process governing the volatility 
will be considered. European call options are then priced under the forward measure as 

SiF,cT,K,t) ^ pit,T)E{{FT - K)+\Ft ^ F,at ^ a), (2) 

where p(t, T) is the price at t of a zero-coupon bond paying 1 at maturity T (see [2j). 



2 Volatility Swaps 

As it became obvious that volatilities and variances for various underlyings behave in a stochastic manner 
on the market, several derivative instruments dealing with those quantities where introduced. 



2.1 Principle of The swaps 

A variance swap consists of a strike K and a realized average variance 1/T J'°+^ afdt over some time- 
period to, to + T. The pay-off function is then the difference of these two quantities. There are known 
closed form solution for variance swaps, but it seems, that these instruments are rather volatile and risky. 

A volatility swap is a bet on a realized averag ed volatility J^°^^ a^dt over some time-period + 

T. By taking the square root of the realized variance, this instrument has a somewhat less volatile behavior 
than a variance swap. But mathematically, it is exactly this square root that makes computations more 
difficult. 

Here we propose to compute exact solutions to the fair value of a volatility swap 



under the SABR model with t E [to,tQ + T]. In the following, we shall write 




a^dt\at = cr), 

so that the fair price of a volatility swap becomes p(t,T){Kt — K). The task is now to compute kj. 
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2.2 Getting Rid of The Square Root 

It is a well-known fact, that for y S K"*" 



^ ' r '-«p(-»'ri ^ (3) 



we therefore obtain 



1 / /■*°+^ 



1 ^ p l-exp(-x^/;;+^a,^d^) , 

■(/ 12 kt=«^) 



T0F 







We will hence have to compute the quantity 



1 pi- E(exp(-a;2 a^d^lf^f = 

io a;2 



(4) 



/■to+T 

/(t, a, x) := E( exp(-a;2 J a'^dt)\at = a) (5) 
and integrate along x as in the previous equality. 

Note, that x is not stochastic but a variable of integration. Also, we did not fully justify the interchange 
of the order of integration, which shall be justified a posteriori. Also, if to <t <to +T, then / will also 
depend on the realized volatility Vt ■= cr^dt "up to f". We suppose at first, that to <t <to+T and 
write 



CF^dt)\cFt = cr) 

t-ta+T 

= E(exp(— a;^ft)exp(— a;^ / afdt)\at = (T) 

rto+T 

= exp(— .'j;^z/i)E( cxp(— / afdt)\at — a) 



= exp{-x ut)ip(t,a,x). (6) 

2.3 A differential Equation for cp 

The Fenman-Kac formula gives the correspondence between expectation computations for a stochastic 
process and differential equations. For a stochastic process governed by the equation 

dXt = ii{Xt, t)dt + a{Xt, t)dWt, (7) 

the function 

u{x,t):=E{J^ eM- ViXr)dT)f{X„s)ds + exp{- V{X,)ds)-^{XT)) (8) 

is the solution to the differential equation 

(7^ fa; t) 

dtu{x, t) + ijl{x, t)d^u{x, t) + y dlu{x, t) - V{x, t)u{x, t) + f{x, t) = 0, (9) 
subject to the terminal condition 

u{x,T) = '^{x). (10) 
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In our situation, we have fj,{at,t) = = f{(Jt,t), 'ii{x) — 1, V{at) = x^af and the stochastic diffusion 
term is a<7t, so that ip{t,a,x) is, for a given x, solution to the differential equation 

2 2 

dtip{t, a, x) + ^^^dlip{t, cr, x) - x^a^ifit, cr,x)^0 (11) 

2 2x^a'^ 
<^ -—dt'pit,(T,x) = a'^dlLp{t,a,x) —Lp{t,a,x), (12) 

subject to the terminal condition 

^(T,a,x) = l. (13) 

Note, that here, x is only a parameter, over which we will have to integrate later on. Here, x is 
indexing a family of differential equations. We shall hence drop it from the variables in ip. Writing 
y := \f2xa I a and ^(t,y) :— ip{t, a{y)), we are to solve 

-^dMt,y) = y^dli^{t,y) - y^i^{t,y), (14) 
subject to the terminal condition 

^(T,y) = l. (15) 
Clearly, we will have (p{t, a, x) = ^{t, ^^^). 

3 Solving the differential equation 

If it wasn't for the missing ydy term, the differential operator appearing on the right-hand side of equation 
[14]would be the modified Bessel differential operator, whose solutions are well-known. In order to retrieve 
this equation, we define g{t,y) :— y^^/'^ip{t,y). This function satisfies the equation 



1 / N 2 

-:9{i.v) 

4 a 

subject to the terminal condition 



-9{t, y) - ^dtg{t, y) = y^d^^git, y) + ydyg{t, y) - y^git, y), (16) 



,iT.y) = ±. (17) 

The solutions to the equation 

k'fiy) - y'd'jiy) + ydyfiy) - y^fiy), (18) 

are the modified Bessel functions of the first and second kind Ik{y) and Kk{y). The functions Ik{y) read 

x-^ 1 ^y\2m+k 



m>0 



The terminal condition y ^1"^ can be expressed in terms of these functions as 

1 1 (-ir(2n-l/2)r(,i-l/2) , 

^ ^ ^ S ' ' 

The convergence of this series is pointwise. We may now write the solution to the equation 1161 as 

y) = E (-l)"('"-)/y("-V2) ,^^^_^^^^^) (^^(j, (21) 
ri>0 ^2".- 

where E„ := ^((2?! - 1/2)^ - 1/4) = a'^n{2n - 1). Consequently, 

i^it, y) = E (-l)"(^"-V2)r(n-l/2) (|^V2,^^^_^^^(^) (^^^(^ ^ _ (22) 

n>0 
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4 Recovering The Square Root 



We shall now compute Kt for a volatility swap under the SABR model, by integrating over x. At first, 
we will suppose to <t < to+T. The case t <tQ will be treated afterwards. We thus have 



^ ^ , 1 - ex.p{-x'^i't)y^{t, a, x) 



TV^ Jo 

1 f°° 1- exp(-a;^ft) + exp(-a;^i/f ) (l - tp{t, a, x)) 



Ty^ Jo 



dx 



If J_ r exp(-a;'i^t)(l - ^{t,cy,x)) ^ 

Ir^, 1 /-exp(-a;2;.,)(l-V;(t,^)) 
If ^, /•~exp(-^j/2)(i_^(^^y)) 



da; 



T L aVTT Jo y J 



To lighten up the notations, we shall write J := ^^p( 2<t^ ^^K^ ^'"^'^^^ dy. Moreover, we write 
_ (-i)"(2n-i/2)r(n-i/2) ^^^^ ^^^^^ _ (|) '^'/2n-i/2 (y)- We wiU Split the integral as 



JO 

Jo + -^oc 



exp(-^)(l - E„>o «ni^n(y) exp(i;„(T + to- t))) 



2 ^2/ 



where 



Jo y 



exp(-^)( En>i«n-^ n(y) exp(^„(T + to - *))) 

V 2m+2n 



Joe = - = -2 dy. 

Jo y 



Observe, that i^„(y) := Em>o m!r(m+2n+i/2) (2) • compute Jq, we insert oq = -l/2r(-l/2) 



0F, and Fo = Em>o nz!r(m+i/2) ( f)^"- Therefore, 



exnC-^^V -1(1/2) fy\'2m 

2—^ — dy 

r ( y%v- -r(l/2) ,y,2„a,y,-2. 



m> 



'^m!r(m + l/2)^2^ 4^2^ 
^ -r(l/2) , 2^-2, 



m> 



Mm!r(m + 1/2) io -2 ^2 



where the interchange of the integral and the sum is justified by the fact, that all terms which are to 
integrate have the same sign, so that the convergence of the integral of the sum is equivalent to the 
convergence of the sum of the integrals. 
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This integral is fairly easily computed: 

1,2 



POO 

= J exp(-,s)(i)"-^/%-3/2d. 



= (|)'"-^/^r(m-l/2). 



We arrive at 



m>l 



4to!(to- 1/2)m' 

where erfi() is the imaginary error function. Hence, we get 

^t = ^{l- ^^F[^ erfi (^) + 1 - exp(./4)] + ^J^} 
We now procede to compute the other terms in Joo- There, we need to compute for n > 1 

Jo y' ^ ~ Jo z> ^mir{m + 2n + 1/2)^2' 



m>0 ^ ' ' ^ 

,2 



El f y N/y\2m+2n-2 

, 4m!r(m + 2n+ 1/2) 70 ""^(-^^(2) 



^ 3^ f°° / y \jy\m+n-3/2ydy 

^^4mT{m + 2n+l/2) Jo ""^^ z 4 ' 2 

Er(n + m— 1/2) /2sm+n-l/2 
4m!r(m + 2n+ l/2)'4'' 

m>0 ' 

r(n-l/2)(n-l/2)(") .z.m+n-i/2 
^^4TO!r(2n+ l/2)(2n+ l/2)('") H'^ 

- ^("-'/'^■(^)"-^/V,(n-l/2;2n + l/2;./4), 



4r(2n+l/2)M^ 



where iFi{a;b;z) is the confluent hypergeometric function. Inserting the exponential factors and the 
a„'s, we arrive at 

7 \- ^p^T^. /■°" exp(~i^)f„(j/) ^ 

n>l •''^ ^ 



4n!r(2n-l/2) 
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Note, that Jo may also be expressed as a confluent hypergeometric function, as 

° ^ 4m!r(m + l/2) 

-r(l/2) r(m-l/2) .Z.m 

A ^A' ^^m!r(m + l/2)'4' 



\2 



.™(£)-(..,-V.;.V.;j.-i 



We finally arrive at 



^ ^ 6„e^"(^+*»-*)C"if^i (n - 1/2; 2n + 1/2; C), (23) 



n>0 



where 



(-l)"+ir(n- 1/2)2 



^" 20Fn!r(2n- 1/2) ' 

En = a2n(2n-l), 
rto+T 

vt — a^dt is the realized volatility up to t, 



c 



2a^vt 4' 

(7 = is the volatility at time t and i-Fi(; ; ) are confluent hypergeometric functions. 
The fair value of a volatility swap with <t <tQ +T is hence 

P^t^ T) E &ne^"(^+*'-*)C"ii^i(" - 1/2; 2n + 1/2; O-k). (24) 

n>0 

4.1 Value at t = to + 7"- 

We replace t with T + io- The formula [23] simplifies to 



«*o+T = ^{ E ^"C"i^i(" - 1/2; 2n + 1/2; C)}. (25) 



n>0 

Inserting the series expansion of the confluent hypergeometric function, we arrive at 



= V^/ V (-l)"+^r(»-l/2)^ (n-l/2)('») ^,^„,^ 
T I ^„ 2v^n!r(2n- 1/2) (2n + l/2)Mm!^ / 



n,m>0 



\M r v;- (~l)"+i(2n - l/2)r(n + m - 1/2) r(n - 1/2) ^^^^ 
T I ^ 2^/^n!r(2n + m+ 1/2) m! ^ /' 

n,m>0 

By regrouping the terms with same exponents, we get 



V^rV.. r(g-l/2) ^ (-l)"+i(2n-l/2)r(n-l/2) ^ 

to+T T ^ n!(s-7i)!r(s + n+l/2) /' 

s>o ^ o<n<s ^ ; V I I / ; 

The first term, i.e. the constant multiplying ^ is 1. The second, with s = 1 and n = 0, 1 is equal to 

-(-i/2)r(-i/2)/r(3/2) + (3/2)r(i/2)/r(5/2) ^ -r(i/2)/r(3/2) + r(i/2)/r(3/2) = o. 

For s = 2 and n = 0, 1, 2, one obtains -(-l/2)r(-l/2)/2r(5/2) + (3/2)r(l/2)/r(7/2)-(7/2)r(3/2)/2r(9/2) : 
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-r(l/2)/2r(5/2) +3r(l/2)/5r(5/2) - r(3/2)/2r(7/2) = (-5r(l/2) + 6r(l/2) - 2r(3/2))/10r(5/2) = 0, 
and one may show by a recurrence relation, that ah these sums equal zero, we thus arrive at 



which is the expected result. 

This is actually the boundary value for the functional calculus we are going to study in the next section 



5 Verification By Functional Calculus 

In the computation of the solution, we somewhat ignored the mathematical difficulties which could arise 
from the exchange of orders of integration, especially when taking expectations and integrating over a 
parameter was involved. In this paragraph, we investigate the solution (|23p in the light of Rama Cont's 
papers |3] & 011 functional stochastic processes. As given by Kt is certainly well-defined. It can 
even be shown, that it is an analytical function. The question is, if it is truly part of the solution of the 
problem at hand. Observe, that 

rto+T 



I to 

is in fact an integral depending on the path chosen for at- As explained in Rama Cont's paper on 
functional calculus, we must have 

DtKt + ^a^a^iVl)^Kt = 0, (26) 

where Dt — dt + and Vjj is the horizontal derivative, while is the vertical derivative. In our case, 
V^C = 0, J^°'^'^a^ds = a^e{t-to)9{T + ta-t), J^"^'^ a^.ds = 0, ^ Kl^^t and the usual chain 
rules apply. We also should have Kt„+T = \/i%+T, which was already verified. 

Let us lighten up the notations and put t := T + to — t and /n(C) '■= iFi{n — 1/2; 2n + 1/2; Q. Kt 
rewrites as: 

T 



n>0 

Computing derivatives, we get 

DtKt = 5tKt + V^Kt = >^{-^6„£;„e^"^C"/«(0} 

n>0 



n>0 



CfLiO+nfniC)]} 



2a^^t{Y.Cfn{0bn 

7l>0 



En 

2^ 



Cn) 



on the other hand. 



[-^tW: 



c/;(c)+n/n(c)]} 



Tl>0 



^ { E bne^"^C-' [2nC/;(C) + n{n - 1)/„(C) + emo] }• 



ra>0 



Now, use the fact, that the hypergeometric functions are solution to the differential equation 

ZiF{'{a;b;z) ^ [z - b)iF{{a;b; z) + aiFi{a;b; z). 
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In our case, 
This yields 



cm) = C(C - 2n - l/2)/;(C) + C(" - l/2)/„(C). 



(V^^K, = 2a2^{^6„e^"-C"[cV;(C) + (C + n)(n-l/2)/„(C)]}. 



By replacing En with their values Q;^rt(2rt — 1), we see that the functional equation (pS)) is satisfied, which 
justifies the equation ([M)) as the solution to the fair value for volatility swaps. 



6 Conclusion 

The fair value of a volatility swap under the SABR model and with t E [to,to + T] has the solution 
. The relation 

^ (-l)"+i(2n-l/2)r(n-l/2) ^ 
ojl^s n!(s-n)!r(s + n+l/2) 

was not yet proven but can be accepted by verification. As a next computation, one may look for the 
value of a volatility swap for t < Iq. The value for a volatility option is also a question that bears some 
interest and should be addressed in the future. 
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